On a new fixed point of the renormalization group operator for area-preserving maps 
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The breakup of the shearless invariant torus with winding number uj — v2 — 1 is studied nu- 
merically using Greene's residue criterion in the standard nontwist map. The residue behavior and 
parameter scaling at the breakup suggests the existence of a new fixed point of the renormalization 
group operator (RGO) for area-preserving maps. The unstable eigenvalues of the RGO at this fixed 
point and the critical scaling exponents of the torus at breakup are computed. 



Area-preserving nontwist maps are low-dimensional 
models of physical systems whose Hamiltonians locally 
violate a nondegeneracy condition (see below) as de- 
scribed, e.g., in Refs. 0, 0, Some applications are 
the study of magnetic field lines in toroidal plasma 

devices jl H 0,0, H, H,E3 and stellarators[llE3 (plasma 
physics), and traveling waves. [l^ coherent structures, 
self-consistent transport (fluid dynamics), and par- 
ticle accelerators. |l5j Nontwist regions have also been 
shown to appear generically in the phase space of area- 
preserving maps that have a tripling bifurcation of an 
elliptic fixed point. [lfil lr?j] Additional references can be 
found in Refs. 0,0 

Of particular interest from a physics perspective is the 
breakup of invariant tori, consisting of q uasiperiodic or- 
bits with irrational winding number. |3flj that often cor- 
respond to transport barriers in the physical system, i.e., 
their existence determines the long-time stability of the 
system. In nontwist maps, the invariant tori that appear 
to be the most resilient to perturbations are the so-called 
shearless tori, which correspond to local extrema in the 
winding number profile of the map. 

Invariant tori at breakup exhibit scale invariance under 
specific phase space re-scalings, which are observed to be 
universal for certain classes of area-preserving maps. To 
interpret these results, a renormalizationgroup frame- 
work has been developed (see, e.g., Refs. I2T HH HH. Eoh. 
For twist maps, it is well understood which fixed point, 
cycle, or strange attractor of the renormalization group 
operator (RGO) is encountered within a given class of 
maps, depending on properties of the winding number of 
the critical torus (see Ref. for a recent review) . For 
nontwist maps, however, only results for the single class 
of shearless critical noble tori, i.e., shearless critical tori 
with winding numbers that have a continued fraction ex- 
pansion tail of l's, are known. The result reported in 
this letter represents the first new fixed point for non- 
twist maps. 

A tool for studying the breakup of a torus with given 
winding number is Greene's residue criterion, originally 
introduced in the context of twist maps.ji^] This method 
is based on the numerical observation that the breakup 



of an invariant torus with irrational winding number uj 
is determined by the stability of nearby periodic orbits. 
Some aspects of this criterion have been proved for non- 
twist maps. [23| 

To study the breakup, one considers a sequence of pe- 
riodic orbits with winding numbers q n /p n converging to 
uj, limn-joo q„/p n = The elements of the sequence 
converging the fastest are the convergents of the con- 
tinued fraction expansion of uj, i.e., [n] '■— q n /p n = 
[do, ffii, • • • , aj, where 



uj = [ao, ax,a 2 , ■ 
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The stability of the corresponding orbits is determined 
by their residues, R n = [2 - Tr(DM Pn )]/4, where Tr 
is the trace and DM Pn is the linearization of the p n 
times iterated map M about the periodic orbit: An or- 
bit is elliptic for < R n < 1, parabolic for R n = 
and R n — 1, and hyperbolic otherwise. The conver- 
gence or divergence of the residue sequence associated 
with the chosen periodic orbit sequence then determines 
whether the torus exists or not, respectively: If the tu- 
torus exists, lim n ^oo \R n \ = 0; if the w-torus is destroyed, 
linin^oo \R n \ = 00. At the breakup, different scenarios 
can be encountered, depending on the class of maps and 
winding number of the invariant torus considered. 

In nontwist maps, the residue criterion was first 
used in Ref. to study the breakup of the shearless 
torus of inverse golden mean 1/7 = (v5 — l)/2 = 
[0, 1, 1, 1, . . .] winding number in the standard nontwist 
map. The residue sequence was discovered to con- 
verge to a six-cycle. Similar studies were conducted for 
other noble shearless tori of winding numbers uj = I/7 2 
(Refs. 0), ^ = [0,2,20,1,1,...] (Ref. EjJ, and 
uj = [0, 1,11, 1,1,1,...] (Ref. 2(3), and the same six-cycle 
was found for all tori in the same symmetry class. 

In this letter we study the breakup of the 
uj = s/2 — 1 = [0, 2, 2, 2,2,.. .] shearless torus, which is 
an example of a non-noble winding number and leads to 
the discovery of a new fixed point of the renormalization 
group operator of area-preserving maps. In contrast 



2 



to twist maps,|27j the periodicity of the elements of 
the continued fraction expansion has not been linked 
to the periodicity of the critical residue sequence and 
therefore our new result would not have been predicted. 
The numerical methods we use and their accuracy are 
discussed in Refs. 0, S H|| and we refer the reader to 
these publications for details. 

As our specific model we use the standard nontwist 
map (SNM) M as introduced in Ref. El 

x l+ i = x l + a(l- yf +1 ) 

Vi+i = Ui - 6sin(27TXj) , (2) 

where (x,y) £ T x I arc phase space coordinates and 
a, b £ R arc parameters. This map is area-preserving and 
violates the twist condition, dxi + i(xi,yi)/dyi ^ 0, along 
a curve in phase space. Although the SNM is not generic 
due to its symmetries, it models the essential features of 
nontwist systems with a local, approximately quadratic 
extremum of the winding number profile. 

One important characteristic of nontwist maps is the 
existence of multiple orbit chains (up and down orbits) 
of the same winding number, which can undergo bifur- 
cations when the map parameters a and b are changed. 
When two invariant tori collide, the winding number pro- 
file shows a local extremum and the orbit at collision is 
referred to as the shearless torus. For a given winding 
number these collisions occur along bifurcation curves 
b u (a) in parameter space. Q 

In order to study a shearless invariant torus, its 
bifurcation curve is found numerically by approximating 
it by the bifurcation curves, 6r n i(a), of nearby periodic 
orbits with winding numbers that are the continued 
fraction convergents of to. Greene's residue criterion 
can then be used to determine where on b u (a) the 
shearless torus still exists: At parameter values a and 
the best known approximation to &tj(a), the residues of 
all periodic orbits of convergents that have not collided, 
here the orbits [n] with even n, are computed. Their 
limiting behavior for n — > oo reveals the status of the 
torus. By repeating the procedure for various values 
of a, with alternating residue convergence to and oo, 
the parameter values of the shearless torus breakup, 
(a c ,b u (a c )), can be determined to high precision. 

The study of the breakup of tori with non-noble wind- 
ing numbers is difficult because, due to numerical limita- 
tions, only periodic orbits for a small number of elements 
of the continued fraction expansion can be found. There- 
fore one cannot numerically distinguish between a torus 
that has 2's as all entries in the continued fraction expan- 
sion and one that has 2's until one reaches the numerical 
limit, and then l's for the tail (i.e., a noble number). 

To make a definite prediction for the oj = v2 — 1 
torus we study the breakup of a series of sixteen in- 
variant tori Ti, starting with u> = [0,1,1,1,...] up to 
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FIG. 1: 2-cycle of a c ,i and b c ,i differences in approximating 
the critical shearless torus lj = y2 — 1 by noble tori Ti . The 
average (negative) slope has been added to the data. 



uj = [0, 2, . . . , 2, 1, . . .] with fifteen 2's in the continued 
fraction expansion. 

The breakup parameters (a c ,j> ^c,i) for the tori Tj 
are found to converge exponentially to the critical 
parameters (a c ,oo> b C:00 ) of the lu = \/2 — 1 torus. 
Plotting the logarithms of the differences |a C) » — a c ,oo| 
and \b Cl i — 6 c ,oo| ; corrected for their average slopes of 
c a = -0.8986 ± 0.0031 and c b = -0.8984 ± 0.0002, 
respectively (see Fig. 0, one observes a period-two 
oscillation. This result should be compared to Fig. 5 
in Ref. I^, where a similar study was conducted for 
the (one-parameter) standard twist map. The plot of 
parameter vs. i showed a straight line with negative 
slope. 

Figure [21 depicts the behavior of a few of the criti- 
cal residues (residues at breakup) on the s± symmetry 
line[3l| as a function of [n]. The values for different 
tori are shifted by 4 along the y-axis to avoid overlap. 
The results show that the more 2's are included, the fur- 
ther the familiar six-cycle for noble winding numbers gets 
pushed towards higher n values. We conjecture that the 
emerging pattern for small n values represents the crit- 
ical residue pattern of the u> = [0, 2,2,.. .] torus. The 
critical residues at 

(a c ,ooA,oo) = (0.446710414656,0.838135537624831489) 

along si converge to the single value R c = 0.621723 for 
the down orbits, and the single value R c = —0.909118 for 
the up orbits. The same result is found along the other 
symmetry lines. 

As discussed in Refs.QIIH close to the critical breakup 
value, the 6r„i (a) obey a scaling law 

&[„] = b c ,oo + B{n) <5r" , (3) 

where B(n) is numerically found to be periodic in n as 
n — > oo with period 4. This period is the same as the 
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FIG. 2: Residue behavior at the shearless T, torus breakup 
for down orbits on s\. For clarity the residues of the tori have 
been shifted upward by 4. (See text for definition of Too.) 



period of the critical fixed cycle of the renormalization 
group operator 1Z (RGO), i.e., a critical fixed point of 

ft 4 , on 

The scaling can be observed by plotting 



In (b [n 



+i] 



B(n) — n In Si , 



where B(n) = ln(B(n + l)/<5i — B(n)) is also peri- 
odic in n. This is shown (on si) in Fig. [21 where 
for clarity only the offsets of In — &[„]) about 

the average slope are shown. The slope was calcu- 
lated from the last 16 difference values by averaging the 
last 12 slopes [ln(6[„ +5] - b [n+4] ) - ln(6 [n+1] - b [n] )] /12, 
with n = 6, . . . , 17. The periodicity of B(n) makes it 
possible to obtain a better approximation for & CjOC (see 
Ref.Q) 



b c ,oc ~ b[22] 



(b[22] - b[i8]) (^[22] - ^[21]) 
(k[19j _ & [!8]) ~~ ( & [ 22 ] ~~ ^[21]) 



This value is used to find the critical residue pattern la- 
beled Too in Fig. Compared to the critical residue 
six-cycle for noble winding numbers, the four-cycle of 
6r„i differences indicates that the to = \[2 — 1 breakup 
exhibits a residue two-cycle along each symmetry line 
(only half of the periodic orbits exist at breakup). As 
discussed above, in each two cycle, the numerical values 
of the residues are found to be the same. 

As in previous studies, the critical torus exhibits in- 
variance under local re-scaling of the phase space in 
the neighborhood of the symmetry lines. Following 
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FIG. 3: 4-cycle of 6[ n ] differences in approximating the critical 
shearless uj — \p2 — 1 torus at a c ^ = 0.446710414656 on si. 



Refs.liHili we compute the scaling factors a and (3 such 
that the torus in the vicinity of its intersection with S3 is 
invariant under (x',y') — > (a x', (3 N y'). (The exponent 
N is the length of the critical cycle, N — 4 for uo — — 1 
and N = 12 for noble winding numbers.) These factors 
are found from the limiting behavior of convergent peri- 
odic orbits: Denoting by (x' n ±,y' n ±) the symmetry line 
coordinates jlgj of the point on the up (+) or down (— ) 
orbit of the [n]th convergent that is located closest to 
(0,0), we compute (see Tabled 
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a± 



lim 



l n,± 



b n+N,± 



, p% = lim 



y' n ,± 



Vn+N,± 



(4) 



The scaling invariance of the torus at breakup and of the 
nearby periodic orbits is illustrated in Fig. ^ 

As shown in Ref. Ilfll the numerical data obtained can 
be used to compute the unstable eigenvalues, 5i and 62, 
of the RGO K by 



and 



Tn = lim 

0^ n— >oo 



TW = lim 

On n— »oo 



-TV] («c) 



b [n] (ac) - b c 



[n+N] 
*c [n] 



(5) 



(G) 



where a c [„] is the a value at which the w-torus breaks up 
along the &r n ](o) bifurcation curve. 

Our results are displayed in Table For comparison 
we also show the corresponding values for the breakup of 
tori with noble numbers. 

In summary, we found a new critical cycle of the renor- 
malization group operator 1Z, a fixed point of 1Z 4 , which 
governs the breakup of the ui = y/2 — 1 shearless invari- 
ant torus. In analogy with the breakup of tori with noble 
winding numbers, we expect this result to be the same 
for all shearless tori (in the same universality class) with 
continued fraction expansion tails consisting of 2's. 
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[0,2,2,2,2,2,...] 

dawn [7] + 
up [7] □ 



-0.012 - / 
-0.014 



-0.015 -0.01 -0.005 



0.005 0.01 0.015 



0~$La [0,2,2,2,2,2,...] 
.+ +. down [1 1] 



up [1 1 ] □ 



-0.0004 -0.0003 -0.0002 -0.0001 0.0000 0.0001 0.0002 0.0003 0.0004 



FIG. 4: Two levels of magnification in symmetry line coor- 
dinates (x',y') of the uj — \/2 — 1 torus at breakup. Also 
shown are the nearby up and down orbits of the [7]th (top) 
and [ll]th (bottom) continued fraction convergents. 



tail 


[...,1,1,1,...] 


[...,2,2,2,...] 


cycle 


12 


4 


a 


1.6179 


2.4725 


(3 


1.6579 


2.8146 


Si 


2.680 


6.311 


52 


1.584 


2.455 



TABLE I: Universal breakup values for noble tori (N = 12) 
and uj = \p2 — 1 (N = 4) . Values for noble tori are from 
Ref.EE 
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